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Abstract 
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• Programming language used: GNU MAKE, DiANA, FoRM, FORTRAN?? 

• Memory required to execute with typical data: Up to about 10 MB 

• No. of processors used: 1 

• No. of bytes in distributed program, including test data, etc.: 369516 bytes, 
including tutorial in postscript format 

• Distribution format: tar gzip file 

• High-speed storage required: from 1.5 to 30 MB, depending on modules 
present and unfolding of examples 

• No. of cards in combined program and test deck: 

• Keywords: Radiative corrections, automatic calculation, Feynman diagram, 
computer algebra, AlTALC, DiANA 

• Nature of the physical problem: Calculation of differential cross sections for 
e"^e~ annhilation in one-loop approximation. 

• Method of solution: Generation and perturbative analysis of Feynman dia- 
grams with later evaluation of matrix elements and form factors. 

• Restriction of the complexity of the problem: The limit of application is, 
for the moment, the 2^2 particles reactions in the electroweak standard 
model. 

• Typical running time: Few minutes, being highly depending on the complex- 
ity of the process and the Fortran compiler. 



2 Introduction 

aITALC is designed as a tool to perform automated perturbative higher order 
calculations of cross sections in high energy physics. 

The main goal is to create numerical programs directly from Feynman rules. 
It is actually developed for few known models like the electro-weak standard 
model (EWSM) and quantum electro-dynamic (QED), and the hmit of ap- 
plication is, for the moment, the 2 — > 2 particles reactions involving only 
external fermions including the one-loop virtual corrections and the soft pho- 
ton bremsstrahlung. 

It is expected from the user to cook his or her own process within a set of basic 

and/or advanced ingredients. As a reward to such effort comes the OTitput of 
aITALC: Fast and reliable FORTRAN code with differential and integrated 
cross sections, analytical expressions for the transition amplitude and nice 
drawings of the individual Feynman graphs. 

The tool was intended to integrate only free-of-charge packages existing on the 
market. As the calculation proceeds, in a modular fashion, it profits from the 
Diana package [1] (based itself on the Qgraf [2] code) for the generation and 
analysis of Feynman graphs, the Form [3] language when dealing analytically 
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with the large expressions in the amphtudes and, finally, the LoopTools [4] 
library (also integrating the FF package [5]) for the numerical calculation of 
the loop integrals. 

This package is developed from the evolution of the static Fortran code 
TOPFIT [6], intended for the calculation of top-pair production. Some of the 
features then available such as the hard-photon bremsstrahlung or separate 
weak-QED contributions were discarded in pro of a full automation, higher 
technical precision and fiexibility for many other processes. 



3 Technical requirements and configuration 

We present a list of the necessary conditions to install and run aITALC on 
a computer. 

System compulsory requirements: 

■ A computer under Linux ^ operating system with standard ^ 

GNU Make utility 
C compiler 

Fortran?? compiler (with preprocessing ability) 
Suggested optional features: 

■ CPU speed > 333 MHz 

■ RAM Memory > 128 MB 

■ At least from 1.5 to 10 MB ^ free space on disk for installation 

■ Further limitations given by the individual modules 

The package is licensed under the GNU general public license and can be found 
under the official web page |http : //www-zeuthen . desy . de/theory/aitalc 

The documentation [?] can be accessed within the package distribution in 
postscript format or at the web page also as pdf . 

The installation includes an accommodation to your system settings via con- 
f igure and the later extraction of individual packages (if needed) and exam- 
ples. This process is also controlled by a make instruction, conforming to the 
standards of the Unix packages. 

^ InstaUation on other Unix systems might be possible if there are aheady some 
modules present. The package has been successfully build in Solaris for testing 
purposes. 

^ Free packages can be obtained from 
and gcc. 

^ The fewer modules you have already installed the larger the size for including 
them. Another extra 30 MB will be needed when building the three examples pro- 
posed. 



http : //www . gnu . org under the names make 
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Fig. 1. Logical flow chart structure for AlTALC. 
4 Makefile like an environmental language 

Automation requires effective inter-communication and a smart way to or- 
ganize the tasks. This second concept has been widely exploited through the 
Makefile environment, being its conceptual organization a pro in order to build 
different sections of the calculation and provide a chain dependency between 
files and modules. 

The Makefile environment accomplishes 

• The simplification of user interface 

• Building the sections tree, loop and f ortran 

• Running the driver process file 

• Writing the intermediate information Diana Form — > Fortran in a 
row 

• Compiling the numerical program leading to the final results 

in a modular fashion. The execution is carried out sequentially following the 
control flow structure designed in Fig. 1. 

Despite of the simplicity, the idea of calculating a process by a make instruction 
has another major advantage: reproducibility, being understood as the abil- 
ity to recreate once and again similar processes for testing or cross-checking 
purposes. On the contrary, customizing of an specific process turns out to be 
a task spread among different parts: driver file, add-on files for neglecting the 
default ones and specific commands in the main fortran programs. 

We will dedicate next parts to a deeper discussion of how each module ad- 
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dresses its tasks and the description of the modus operandi. 



5 Feynman diagram generation with Diana 

Let us introduce some basic concepts about DiANA. If we are thinking of some 
process, a very intuitive approach is to draw some Feynman diagrams to get 
a feehng and support the calculation which comes later. 

Such a task is very likely to be implemented with the help of computers since 
it is basically a combinatorial problem. Once we define a model, the external 
particles, and the level of loops, the problem can be solved. Indeed, DiANA 
will do it for us. 

5.1 Basic driver file: Input 

The driver file (e.g. process . ini) serves as basic input. It contains primarily: 

• physical model, 

• incoming and outgoing particles, 

• process control flags and options. 

These items are used in the Diana "create" file, that will serve as main user 
input stream for DiANA at each level, i.e, tree or loop. Only the number of 
loops is left to be arranged and written by the Makefile on the "create" . For 
example, it is compulsory to check and agree with the five lines 

SET _processnaine=muon_production 

\Begin (model ,EWSM. model) 

ingoing le( ;pl) ,Le ( ;p4) ; 
outgoing lm( ; -p2) ,Lm( ; -p3) ; 

options=notadp , onshell ; 

in process.ini before creating any process. For the case of /x-pair produc- 
tion, as considered in the previous lines, the full driver file is shown in the 
Fig. 2. Only one remark at this point: special care is needed with the mo- 
menta definitions. They all should be assigned to be incoming and clockwise 
indexed (e.g. pi addressed to be the first incoming particle, -p2 to be the first 
outgoing one and so on), so we ensure a matching with the convention used 
in the predefined topologies. 
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Driver file for Diana 2.35 <process , ini> 
This file is part of alTALC 

Copyright (C) 2004 Alejandro Lorca <ale jandro . lo 
rith your taste, mainly 



DIANA RUNNING SETTINGS 



■ca@desy . de> 
ifter the 



■ MODIFY HERE! 



oduction 

n loop momenta and naked integratit 



'SET _MARk\cOP'"YES" 
'SET rim - "rim" 
" Next (YES/NO) avoid/load 

** generic topology tables and internal abstract momenta 

SET _NO_TRBLES="YES" 

*SET _NO_TABLE_MOMENTA="NO" 

" To be adjusted bv autoconf 
SET _DIANANAME-diana35 

" if yo"^ want to edit the topologies: 
'SET _TC?0LOGYEDITORNAME=tedi 

or.ly interpret 
\openlanguage (create .tml) 
\MINVERSION (2 . 35) 



\Begin (program) 

' MODEL TYPE 



" MODIFY HERE! " 
*\Begin (model, QED .model) 
\Begin (model, EWSM, model) 
\End(model} 



DESCRIPTION OF OUR PROCESS (INGOING AND OUTGOING PARTICLES) 



\Begin (process) 

** MODIFY HERE! ** 

ingoing le ( ;pl) , Le (; p4 ) ; 

outgoing Im (; -p2) , Lm (; -p3) ; 

** To be adjusted in the Makefile 

loops=N; 

SEnd (process) 



process.ini Page 2/2 




Mcopmomenta (q) 












* For the topology tables 






*\ token (bridge, b) 






*\token (chord, k) 






*\token (loop, q) 






\include (runningrestrictions .prg) 






* LOADING/ GREAT ION OF TOPOLOGIES 






" When creating a topology from the 


scratch. 


please save it into 


" a file and rerun Diana with that topology 


file. 


" To be adjusted with the Makefile 






*\tables (41egsNloopS .tab) 






\begin (topology, 41egsNloops . top) 






\End (topology) 






\beqin (agraf opt ions ) 






" MODIFY HERE! " 












opt ions -r.ctadp, onshell ; 






\End(qgrafoptions) 






' POSTSCRIPT FLAGS AND OPTIONS 


\begin (header) 


" MODIFY HERE! " 






SET MakePs - " ! " 






SET Makeinf = " ! " 






SET MakeEps = " ! " 






SET Makeinf oEps - "I" 






debug off 






forced execute 






extra call 






messages disable 






\End (header) 






' TML LANGUAGE 


\ include (thetmlstuff ,prg) 


\ indices (mul, mia2, mu3, mu4,mia5, mu6, mu7 


mu8,mu9. 


mul0,mull,mul2) 


\vectors(p,pl,p2,p3,p4,p5,p6,k,kl,k2 


k3,k4,k5 


,k6, q, ql, q2, zero) 


\End (program) 



Fig. 2. Snapshot for process . ini driver file. 

The comment character is * and the ** MODIFY HERE! ** sentences indicate 
where some modifications arc accepted. Regarding the flags for the initial 
settings, the topology editor (Tedi) can be launched at running time by un- 
commenting the line: 
*SET _TOPOLOGYEDITORNAME ="tedi". 



5.2 Implemented models 



Once the driver file is fully set up, DiANA is able to offer some of the in- 
formation to Qgraf and, by looking to the allowed couplings given in the 
model file, create the possible diagrams appearing at each level. The presently 
available model files are QED. model and EWSM. model. They are summarized 
in Tab. A.l and Tab. A.2. 

Both of these models share basic properties according to their phenomenolog- 
ical content. We can distinguish between the following sections: 

• Propagator definitions, with general structure 

[ Particle , Antiparticle ; Prototype ; Function ( Prop-label , arguments 
)* factor ; Mass-sq. ] 

• Boson, e.g. [Wm,Wp ;W; VV(num,ind:l ,ind:2 ,vec, 2)*i_ ;MMW] 

• Fermion, e.g. [le,Le;l; FF(num,fnum,vec, Mle)*i_; MMle] 

• Ghost, e.g. [ghZ.GhZ; Z ; SS(num,l)*i_; MMZ] 
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• Vertex definitions, as 

[Interacting particles; ; Function (arguments) *factor] 

in which the specific Feynman rules will act later. They are, to some extent, 
part of the model but, from a computational point of view, they belong to the 
module controlled by FoRM programs and therefore encoded in a procedure 
(ctf eynmanrules). 



5.3 Topology editor 



Tedi is designed to allow the user ultimate fiexibility on topology naming, 
shaping and even defining momenta. It was developed together with Diana 
as a tool for dealing with topological issues of the Feynman diagrams, and 
requires the X window server being running in your system. The 'windows' and 
'tab' structure renders an intuitive use, as it is shown in Fig. 3. Furthermore, 
help and printing facilities are also embedded in the program. 



File Edit 



Topology editor: file 4legs1 loops.top 



Controls Help 



EXIT t lO CM ^ 

56 



Shape I Labels | Numbers / Marks | Directions Momenta 



JdU hoo 



61 
62 
63 
64 
65 
66 
67 
68 
69 
70 
71 
IZ 
73 
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[300 !350 Hoo 
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Space to push "Zoom" tool button, Tab/Shifl+Tab - next/prev 



Fig. 3. Snapshot for the momenta tab in Tedi. 
The content of the five tabs has to do with 



Shape: handles with the aesthetics of the Feynman diagrams. Here the users 
can reshape the different legs and loops of the topology according to an 
optimized algorithm or their own taste. Aligning to a grid is also possible 
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by adapting the alignment distance and toggling the key points between 
nodes for each object. 

• Labels: External legs, internal propagators and vertices are given a label in 
Diana. The location of this label might be varied with respect to the graph 
to avoid non-legible superpositions. 

• Numbers/ Marks: Every line has an index that is negative for external legs 
and positive for internal propagators. Working with this tab one can ensure 
to give the highest index values to a loop or to fixed-momenta lines. 

• Directions: The sign of the momentum flow in each internal line can be 
flipped by clicking on the arrow with the right-button mouse. 

• Momenta: Alternative symbols can be arranged for the momenta attribute. 

Apart from this main five features, a bar on the left scrolls between the dif- 
ferent topologies that were generated for the specific process. 

A deeper discussion of the possibilities of Tedi is out of the focus of a standard 
aITALC description. Indeed it has been chosen that the default running of 
any of the examples does not require the topology editor at any moment. We 
will close this summary with the Fig. 4 in order to give the user an impression 
of the standard graphical output of Tedi. 





Fig. 4. Collection of topologies given by Tedi for the Bhabha scattering at one-loop 
level. Tadpoles on vertices account for counterterms diagrams. 
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5-4 Advanced setup 



But what if we want to exclude some couplings or modify the aspect of the 
diagrams? 

After the installation of aITALC it is also possible to use Diana ^ as an 
independent package. We will show shortly the tested modifications that can 
be performed at this level. 

After running an example (e.g. muon_production), inside the tree directory 
(also the loop directory in other examples) it is possible to select some files to 
be considered instead of the default ones coming from the following directory: 
$AITALCHOME/diana/prg/. This is automatically done by placing those files on 
the exainples/muon_production/ directory and remaking the process. When 
the process has to be run for the second time, the instructions are: make clean 
and then make again. 

We do not suggest to copy and modify any other file from prg since they are 
mostly functions using advanced DiANA declarations. The following files are 
welcomed to be modified: 

• part icleaspect . prg: Defines how every propagator and particle label looks 
like in the .ps and . eps files. Syntax is not described (just a coTiple of com- 
ments), but intuitive. The field label CT is assigned to countcrtcrms. 

• runningrestrictions .prg: Fields or couplings can be excluded here, by 
uncommenting or replacing them where desired ^ . 

5.5 What do we get? 

Let's have a look into the level directory (e.g. tree/ or loop/ as in Tab. 1): 

the organizer file is called Makefile. It has the instructions to create, organize 
and fill the different directories by given commands under the MAKE environ- 
ment. This means that in order to build the process, only a make instruction 
is needed. Depending on the process, modifications, hardware and compiler 
the time for each module varies from seconds till several minutes ^ . At the be- 



^ This extends also to FORM and LOOPTOOLS. 

^ Warning: As soon as you may discard some fields or couplings, be aware that 
the renormalization of parameters and fields was fixed and completed within each 
model and does not know about your intentions, so performing a consistency check 

of the results is mandatory. 

^ If large files are to be compiled, occasionally the system might overload. Every- 
thing taking longer than half an hour is suspicious of going awry. 
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Table 1 

Listing for the tree directory inside the muon production example. 

alorcaOlinux: ~AITALCHOME/examples/niuon_production/tree> Is -XFl 
total 183 



H T*T.T V T* — 


VT* — 
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904.8 




1 4 


1 R 


09 




H T'T.T'V T* — 


■VT* — 




9 




th 




rirl" 


1 A 


1 R 

J. >J 


09 


u' u'Tmn /-iTi TiTTiHrn^ti / 


H T'T.TYT' — 








CLJ-Ui- V^CL 




904.8 


rir-i- 


1 4 

X*! 


1 R 

J- o 


09 


Tnf nFPC!/ 








1 


CLJ-UX V^CL 


oil 




Uu u 


X'l 


1 1; 


09 


no r\.c± J.J-C 


— T*T.T— T* — 
J. W J. 


— T* — 
J. 




1 

X 






•^71 




1 4 

XT 




09 


TTin riTi TiT'/^Hn^t'i riTi ^ttF 


— 3rw — 1l~ 


— 1l~ 








th 


90 




14 


X O 


09 


"F oynri "i nn ^nyyon t Qn ?»Tno Q in 

X 6 J. IIIJ. ^Jll^ LLJ. J. Cll U OllCLillC O < J. 11 


-rw-r- 


-r- 




1 


alorca 


th 


3870 


Oct 


14 


15 


02 


muon_production . in 


-rw-r- 


-r- 




1 


alorca 


th 


4531 


Oct 


14 


15 


02 


MaJse_kit . log 


-rw-r- 


-r- 




1 


alorca 


th 


149 


Oct 


14 


15 


02 


Makefile . log 


-rw-r- 


-r- 




1 


alorca 


th 


1088 


Oct 


14 


15 


02 


do_amplitude . log 


-rw-r- 


-r- 




1 


alorca 


th 


245 


Oct 


14 


15 


02 


joinf f . 1 . log 


-rw-r- 


-r- 




1 


alorca 


th 


1059 


Oct 


14 


15 


02 


joinf f . log 


-rw-r- 


-r- 




1 


alorca 


th 


629 


Oct 


14 


15 


02 


j oinkinematics . log 


-rw-r- 


-r- 




1 


alorca 


th 


121102 


Oct 


14 


15 


02 


joinmm. log 


-rw-r- 


-r- 




1 


alorca 


th 


14779 


Oct 


14 


15 


02 


muon_product ion . ps 


-rw-r- 


-r- 




1 


alorca 


th 


18277 


Oct 


14 


15 


02 


muon_productionInf .ps 



ginning and end of each module, time stamps are placed in the Makefile . log 
file at Tab. 1. 

The directory also offers Form processed files that we will discuss later. At 
the moment we can concentrate in just two aspects: 

• Graphical representation 

The global files Sprocessname .ps and $processnameInf o .ps summarize 
all the Feynman graphs appearing at the process. In the EPS/ (and EPS- 
Info/) directories one may find . eps files with all individual diagrams (and 
momenta- label definition). Fig. 5 represents the two types of output for a 
single vertex diagram for the t-channel exchange in Bhabha scattering. 

• Diagram info 

As it comes out of DiANA, the file $processname . in contains the required 
information about both, global process and individual diagrams. Besides 
a compulsory expression for the amphtudc, specially useful for the FORM 
routines are the preprocessed variables defined for each diagram. Below we 
show the connection between the Fig. 5 and the corresponding output in 
Sprocessname . in. 

*— #[ nil: 

****(qgraf number 11) 
**** (diagram 11) 
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a) b) 




Fig. 5. Encapsulated postscripts diagrams for a typical diagram in Bhabha scat- 
tering, a) Simple diagram from EPS/ directory, b) also with momenta distribution 
from Inf oEPS/. 

#define imml "MMle" 

#define ql "-q-kl" 

#define looplinel "1" 

#define iinm2 "MMle" 

#define q2 "-q-k2" 

#define loopline2 "1" 

#define immS "MMA" 

#define q3 "+q" 

#define looplineS "1" 

#define inim4 "MMA" 

#define q4 "+p4+p3" 

#define loopline4 "0" 

#define kl " (+p2) " 

#define k2 "(-pi)" 

#define DIAGRAM "11" 

#define COUNTER "11" 

#define LINE "4" 

#define FERMIONLOOP "0" 

#define TOPOLOGY "4tlc" 

#define CHANNEL "t" 

#define SKELETON "1" 

#define LOOPTYPE "c" 

#define LOOPTYPENUMBER "3" 

#define PROTOTYPE "llaa" 

#define FERMI ONCHANNEL "T" 

#define 11 "1" 

#define 12 "2" 

#define 13 "1" 

#define 14 "2" 

* 
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g [Amplitude, 11] = 

(-l)*F(3,l,mu2,l,0,l)*(-i_)*e*Qle*FF(l,l,-q-kl,Mle)*i_* 
F(4,l,mu4,l,0,l)*(-i_)*e*Qle*FF(2,l,-q-k2,Mle)*i_* 
F(l,l,mul,l,0,l)*(-i_)*e*Qle*F(2,2,mu3,l,0,l)*(-i_)*e*Qle* 
VV(3,mul,mu2,+q,0)*i_*VV(4,mu3,mu4,+p4+p3,0)*i_; 

*--#] nil: 



6 Algebra manipulation: Library in Form 

The form module acts as a bridge connecting the world of the purely symbolic 
representation of the Feynman diagrams, directly read off the lagrangian, and 
the numerical description of a scattering problem of real particles. 

When Diana has finished one of the levels, a complete description of the pro- 
cess is obtained, but still in an encoded language. Different Form programs 
will prepare the amplitudes in order to render numerically evaluable expres- 
sions. Tracing, commuting 7 matrices, applying the Dirac equation to external 
fermions or introducing the Mandelstam variables is part of this module. 

There are basically two kind of programs: The main routines and the proce- 
dures. We describe the four routines: 

• do_amplitude . f rm: Attempts to simplify the individual amplitude for each 
diagram, extracting out the form factors in a given basis of matrix elements. 

• joinkinematics.frm: Defines how the differential cross section should be 
composed in terms of fermion currents (S, T or U), the normalization of the 
incoming fiux in the cross section formula, soft photon emission and the 
definition of the Mandelstam variables. 

• joinmm.frm: The standard set of matrix elements is considered by this 
program, giving, as output, the multiplication of any two elements appearing 
in the amplitude decomposition. 

• joinff .frm: Once all the diagrams were considered, this part pastes to- 
gether the form factor contributions sorted by topologies. The final trans- 
lation into Fortran code is also done at this step. 

6.1 Procedures in kitFOBM.3/ 

The listing of the procedures is long, and the description of each of them is 
shortly presented in Tab. B.l. 

At least we would like to draw the attention to one option regarding the 
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neglect of the external masses of particles (e.g. light fcrmion masses). This 
can be achieved including the adapted Neglectedmasses . inc file into the 
main directory, together with the driver file process . ini. 

Coming back to our previous listing of the directory at Fig. 1, in the files 
do_amplitude . log and joinnun. log we find information about the individual 
evaluation of each diagram and the kinematical factor corresponding to the 
cross product of matrix elements, respectively. Also, extra functions related 
to soft photon emission, flux factors and Mandelstam variables can be found 
in joinkinematics.log. As an intermediate step, individual amplitudes are 
stored in the FF$processaname/ directory. 



7 Getting numbers: the Fortran code 

So far the expressions that were generated contain large amount of analytical 
work on them. Now one has to execute them and evaluate the parameters, loop 
integrals and finally retrieve finite formfactors that, clustered by topologies, 
will allow together with the matrix elements to calculate the cross section for 
a given process. 

The structure of the f ortran module can be outlined as follows: 

• main.f , parameterlist .hf : they both give the user access to modify the 
output structure of the program and the explicit values of the model re- 
spectively. In the main . f file we find many settings as logical flags, input of 
kinematics and output control. 

• kitFORTRAN: consists of the full set of Fortran subroutines and func- 
tions needed to compile the main . f program. We can divide them into two 
categories, according to the role they play during running time. 

• Global are those that were fixed at installation time and considered to 
be process independent. They are stored at $AITALCHOME/f ortran/src 

■ Local are written at running time by the tree and loop modules and 
depend explicitely on the application. They get placed at Sprocessname/- 
f ortran/src. 

• LoopTools: this library is called with the purpose of numerically calcu- 
lating the most of the loop integrals. 

7.1 Understanding and controlling the output 

A sample output are the Tabs. 2 and 3 with differential cross sections (given 
always in pb). The first column corresponds to the beam energy {^/s — 
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setsqrtsman) given in GeV, the second one to the cos^ (setcost), 6 being 
the angle between the three-momenta pi and — P2- Later on we have succes- 
sively the columns of the Born approximation (Born), the interference terms 
of the tree and loop levels (loop) and the soft photonic corrections (soft). 
These three columns (third, fourth and fifth) are summed into the sixth col- 
umn (B+corr) with the finally corrected cross section. The last column loop''2 
contains the one-loop squared terms in the perturbative approach. The max- 
imum energy fraction that the soft photon may gain out of \/s is limited by 
the variable setf racomega. 

Table 2 

File main, log, for the bhabha_QED example. 

#============================================================================================= 

# alTALC Version 1.0 (29.10.2004) by A.Lorca — T.Riemann 

# 

Ssqrtsman cost dcs (Born) dcs (loop) dcs (soft) dcs(B+l+s) dcs(loop"2) 

0.5000000E+03 -.9000000 . 5238736E+00 . 1045142E+00 - . 2193133E+00 . 4090746E+00 . 9235446E-02 
0.5000000E+03 -.5000000 0.6116007E+00 . 1767866E+00 - . 2922801E+00 . 4961073E+00 . 1727483E-01 
0.5000000E+03 0.0000000 0.1172536E+01 . 3974063E+00 - . 6065847E+00 . 9633579E+00 0.4206711E-01 
0.5000000E+03 0.5000000 0.5504406E+01 0.2126643E+01 - . 3064635E+01 0.4566414E+01 0.2405476E+00 
0.5000000E+03 0.9000000 0.1891184E+03 0.8703341E+02 -. 1165000E+03 0. 1596518E+03 . 1089656E+02 
#============================================================================================= 

#Energy given in GeV, Cross sections in pb, setf racomega= 0.1000000 

If the calculation is correct, all the columns should be constant ^ under varia- 
tion of the ultraviolet (UV) parameter (you can check this by turning on the 
fiag luvcheck). The same check for the infrared (IR) behaviour (iricheck) 
will change the explicit values in columns loop and soft, but the sum should 
still be invariant. The last column may also suffer against IR variation since it 
is neither compensated by doiible soft photon emission nor two-loop divergen- 
cies. Those loop~2 numbers are crucial when they become the leading order 
in e.g. fiavour changing neutral current processes (FCNC) as discussed in [8], 
or just as an estimator of the order of magnitude for the error in the next 
level of perturbation theory. When integrated cross sections are under study, 
the integration region is limited by the setlimcost variable, the ics and f ba 
being defined in terms of dcs as 



dcs = - 1 

dcosfc* 



setlimcost 

ics = (7tot 



= j d COS ^ dcs (2) 

—setlimcost 

setlimcost 

f dcs f dcs , , 

/ dcos^ / dcos^ . 3 

J ics J ics 



—setlimcost 



Warning: as a result of numerical variation, the round off of intermediate results 
will modify the last digits in your results. Such variation increases at the collinear 
cases due to cancellations. 
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Table 3 

File main. log, for the leLe.bS example. 

#============================================================================================= 

# alTALC Version 1.0 (29.10.2004) by A.Lorca — T.Riemann 

# 

Ssqrtsman cost dcs (Born) dcs (loop) dcs (soft) dcs(B+l+s) dcs(loop"2) 

0.5000000E+03 -.9000000 . OOOOOOOE+00 O.OOOOOOOE+00 . OOOOOOOE+00 . OOOOOOOE+00 . 65341 18E-08 
0.5000000E+03 -.5000000 O.OOOOOOOE+00 . OOOOOOOE+00 . OOOOOOOE+00 O.OOOOOOOE+00 . 7095421E-08 
0.5000000E+03 0.0000000 O.OOOOOOOE+00 . OOOOOOOE+00 . OOOOOOOE+00 O.OOOOOOOE+00 . 1048690E-07 
0.5000000E+03 0.5000000 O.OOOOOOOE+00 . OOOOOOOE+00 . OOOOOOOE+00 O.OOOOOOOE+00 . 2084942E-07 
0.5000000E+03 0.9000000 O.OOOOOOOE+00 . OOOOOOOE+00 . OOOOOOOE+00 O.OOOOOOOE+00 . 4643956E-06 
#============================================================================================= 

#Energy given in GeV, Cross sections in pb, setf racomega= 0.1000000 

The integration algorithm is based on a Richardson extrapolation to the 
Romberg integration [9] with four steps. The estimated error is supplied in 
short forms in extra columns. If it does not suffice or for cross checking, the 
same subroutine with eight steps (ICS8) can be called instead of the faster 
ICS. A complete survey of variables is presented in Tab. C.l. 

1.2 A greedy possibility: quadruple precision 

Not every compiler^ allows for extended or quadruple precision. This is a 
feature outside of the ANSI Fortran??. Even if aITALC contains some 
code's implementations outside the ANSI syntax, they are quite common (e.g. 
using underscore, or double complex), so decent compilers won't complain. 

But on high level computation, or simply as comparison, switching on the 
quadruple precision could stabilize numerical results or render amazing agree- 
ment with other calculations. 



8 Examples including loop level 

To finish this descriptive chapter, we would like to present the examples avail- 
able with the distribution. 

The package is delivered with three examples that the user may run with- 
out further considerations. In the directory examples/ we have the following 
processes: 

• /i-pair production at tree level: e^e^ A* /^^ 

• Bhabha scattering in QED (see [10]): e~e~^ —>■ e~e~^ 

^ Until now, the use of quadruple precision has been positively tested under Intel 
Fortran compiler (see 

|http : //www. Intel . com/ software/products/ compiler s/flin/noncom. html . 
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• Fermion Flavour Violation example: e e"^ — > 6s 

under the directory names of muon_production, Bhabha_QED and leLe.bS. 
Table 4 

Typical timings for the different modules. The technical specifications for system 
Laptop arc: Intel Ccntrino 1.5GHz cpu, 512MB RAM, Intel ice and if ort version 
8.1 compilers. For Desktop: Intel Pentium III 853MHz cpu, 256MB RAM, GNU gcc 
and g77 version 3.3.3 compilers. 



Module 


muon_production 
Desktop Laptop 


leLe.bS 

Desktop Laptop 


bhabha.QED 
Desktop Laptop 


tree 
loop 
f ortran 


9" 3" 
15" 4" 


4" 1" 
2'20" 52" 
4:35'56" 33" 


30" 9" 
37" 12" 
26" 8" 


Total 


24" 7" 


4:38'20" 1'26" 


1'33" 29" 



Running time expected for the execution of all the examples were studied on 
two different computers. The results are shown in Tab. 4. One can see that 
when the size of the Fortran code begins to be large (i.e. more than ~200kB 
for a single subroutine to be compiled), then the required time for the GNU 
compiler gets too long, compared to the other processes or compiler. 



A Implemented models 
Table A.l 

Particle content of the QED model (QED .model) 



Code Arguments Particle 


Name 


Lorentz type 


Lcpton 






le, Le (;pi) e", e+ 


electron, positron 


fermion 


Gauge boson 






A (fJ-ilPi) 7 


photon 


boson 
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Table A.2 



Particle content of the electro- weak sector in the standard model (EWSM. model) 



Code Arguments 


Particle 


Name 


Lorentz type 


Leptons 










ne , Ne 






neutrino-e, antineut-e 




run , Nm 






neutrino-//, antineut-/i 


fermion 


nt , Nt 






neutrino-r, antineut-r 




le , Le 




e~, e'^ 


electron, positron 




Im , Lt 




fi-, ^+ 


muon, antimuon 


fermion 


It , Lt 






tau, antitau 




Quarks 










u , U 




u, u 


up, antiup 




c , C 


(.;Pi) 


c, c 


charm, anticharm 


fermion 


t , T 




t, t 


top, antitop 




d , D 




d, d 


down, antidown 




s , S 




s, s 


strange, antistrange 


fermion 


b , B 




b, b 


bottom, antibottom 




Gauge bosons 








A 




7 


photon 




Z 




Z 


Z-boson 


vector 


Wm , Wp 




W-, w+ 


W-boson 




Higgs sector 










H 




H 


Higgs 




GO 


(;Pi) 


X 


would-be 


scalar 


Gm , Gp 






Goldstone bosons 




Faddeev-Pop 


ov ghosts 








ghA, GhA 




77-y, rj^ 


photon ghosts 


scalar 


ghZ, GhZ 


(.;Pi) 


Vz, fjz 


Z-ghosts 


with 


ghm, Ghm 




vw-,fiw- 


W~ -ghosts 


fermion 


ghp, Ghp 




r]w+, 'nw+ 


W"'"-ghosts 


statistics 
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B List of Form procedures 

Table B.l. Different procedures used. 



Name Args. 


Description 


LTtoPV 


Translates the indexing of loop integrals 

il Ulli 1 J- J- 1 llltU 1 -L^ 1 CUllVclltlUll. 


aitalcnotation N 


Translates the Form code into Fortran 

leaLiauic. ±iie aig,u.iiieiii iv — i aiiowb loi 
complex masses while A*" = does not 


analyzeterms 


Part of simplif ygrams 


argsymmetries 


Applies symmetry properties on argu- 
ments for loop functions and Gram de- 

tCi llillidiito 


^ CLLl ^ C7 X 11 o 


j.\J L/CllH-.d l^ljlll L/Uol LlUll KJi- 

den(a,b)*(a-b)=l 


CiHoiioxin 


/\ppiieb onibiioim iQcniiLies oy conveii- 
ing all chains of 7 matrices that are con- 
tracted by their dimensional value 


contgammas 


Simplifies contractions in products of 

gamma matrices 


contractepsilon 


It contracts the product of two e tensors 


ct f eynmanrul e s 


Insertion of Feynman rules with countert- 
erms 


def inegrams 


Defines a FORM global variable with the 
explicit expression of a Gram determinant 


derivateself 


Derivates self-energy functions 


dimens i onf our 


To reach the limit D — > 4 taking care of 
the UV-behaviour of scalar integrals 


diracequation 


Applies Dirac equation to available 
spinors 


continued on next page 
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continued from previous page 


Name 


Args. 


Description 


dummy tovar 


N 


Part of simplif ygrams. stands for the 
amount of arguments in dummyN function 


equal index 




Gives the same index to potentially dif- 
ferent contractions keeping them ordered 
by number 


ext ernalmoment a 




To substitute the internal momenta in 
loops with external ones 


gammaStol 




Introduces an e tensor and 75 to remove 
the product of three 7's 


gammaalgebra 




Prepares 7 chain structures for 
diracequation 


gammamunu 




Introduces indices in the vector structures 


gramback 


N 


Part of simplif ygrams. Returns the 
Gram determinant of order A^ if no sim- 
plification was found 


gramsubstitution 


N,w 


It substitutes the w-e\emcnt to match 
posssible order A^ Gram determinants 


identif yintegrals 




Picks up the different loop integrals and 
provides a set of Form $-variables with 
the right substitutions, inserting the func- 
tion LoopIntegral=LI (n) 


integrationde93 




Procedure to integrate 4,3,2, 1-point func- 
tions with internal momenta decomposi- 
tion like in [11] 


invgammamunu 




Undoes the gammamunu effects 


keepUVloops 




Keeps only integrals actually UV diver- 
gent 


loopsymmetries 




Minimize loop integrals in boxes by sym- 
metrizing arguments 


lorentzinvariants 




Adapts scalar product of vectors to 
Lorentz invariant squared masses 


continued on next page 
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continued from previous page 


Name 


Args. 


Description 


masse onvent ion 




To use single variables for the masses with 
M instead of MM outside arguments 


massivef of a 


J, c, N 


* — 111 1 lie C 1 

L-alculates and saves the tormtactors as 
[13]. j stands for topology, c for left-right 
or I-75 and N=1(0) do (not) print the 
formfactors 


movepslash 




Places the desired to the side of fermion 
chain to apply diracequation 


neglectmass 


zero 


Neglects positive powers of terms as indi- 
cated in the file Neglectedmasses . inc. 
zero is the substituted value 


nosymboliccouplings 




To put explicit values of charges and weak 
couplings 


onshell 




Substitutes external momenta for exter- 
nal particle masses and Mandelstam vari- 
ables 


pslashaway 




Reorder away the non desired position of 


pushgammaS 




Pushes 75 to the right in fermion chains 


pushomegas 




Orders and pushes Left or Right projec- 
tors {ojl,r) to the right in fermion chains 


recoverargumentsde93 




Catchs back the arguments of 
loop integrals in combination with 
integrationde93 


reductionDD02 


N 


Reduction algorithm from scalar tensor 
integrals to master integrals given in [14] . 
N = 1 sorts after each integral type, 
iV = does not sort at all 


reductionLT 




reductionDD02(l) with LOOPTOOLS 
notation 


continued on next page 
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continued from previous page 


Name 


Args. 


Description 


simplifygrams 




Tries to remove inverse Gram determi- 
nants by looking to the numerator (time 
and memory consuming!) 


simplihelp 


N,r 


It helps simplifygrams to run sequen- 
tially r times from the highest power of 
one variable for N order Gram determi- 
nants 


storeself 


j, c, 


Stores self energies, j stands for topol- 
ogy, c for {ujl.ojr) or (1,75) and N=1(0) 
do (not) print the form factors 


threegammastoepsilon 




Introduces and removes the e ten- 
sor simplifying chains of 3 7 matri- 
ces that naturally don't dissappear after 
diracequation 


tracef ermiloops 




Traces possible fermion loops in self- 
energies 


transvorlongit 


X 


Sets the longitudinal and transverse part 
of self energies, keepmg orr-shell but with 
mass^ argument x 


unityCKM 




Makes CKMij = 5ij 


usegammaS 




Uses 75, instead of 0;^,^^^ projectors 


useomegas 




The opposite to usegammaS 


vartodummy 




Part of gramback, opposite as 
dummy tovar 
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C Settings for main.f 

Table C.l. Kinematical and variables and logical flags for different settings in 
main . f . 



Variable 


Type 


Def. value 


Meaning 


Kinematical 








icostloop 


integer 


5 


Dimension of setcostarray 


isqrtsloop 


integer 


5 


Dimension of set sqrt smanarray 


setcostarrav 


double 


data // 


Set of default values for the 
setcost (cos^^) do loop in the 
dcs evaluation 


s et sqrt smanarr ay 


double 


data // 


Set of default values for the 
setsqrtsman (a/s) do loop in the 
ICS evaluation 


setf racomega 


double 


O.ldO 


Limit for the maximum soft pho- 
ton energy 


setlimcost 


double 


.9999d0 


Limit of integration for ics and 
fba 


Process Bags 








Irenorm 


logical 


. true . 


Call the renormalization subrou- 
tine to calculate the counterterms 

parameters 


Iwidth 


logical 


.false . 


Call the loop integrals with com- 
plex values for boson masses. 
For the four-point integral, only 
the case with one massless boson 
(photon) is considered. 


lidentCKM 


logical 


. false . 


Brings the CKM mixing matrix 
in a diagonal form, so no mixing 
occurs 


continued on next page 
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continued from 


previous page 




Variable 


Type 


Def. value 


Meaning 


Checking flags 








luvcheck 


logical 


. false . 


Perform a shift in the mudim 
and delta parameters in LOOP- 
TOOLS 


lircheck 


logical 


. false . 


Jr^ertorm a shift m the lambda pa- 
rameter in LoopTools 


Output flags 








Icostloop 


logical 


.true . 


Performs a do-loop over setcost 


Iprintics 


logical 


.false . 


Performs a do-loop over 
setsqrtsman and prints the 
results of ics 


Iprintf ba 


logical 


. false . 


Performs a do-loop over 
setsqrtsman and prints the 
results of f ba 


llongoutput 


logical 


. false . 


Prints the cross sections in double 
precision format 
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